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Abstract 

We present the method of calculation of the radiative decays of composite quark- 
antiquark QQ systems, with different J PC : (QQ)i n — > ^{QQ)out- The method is 
relativistic invariant, it is based on the double dispersion relation integrals over the 
masses of the composite mesons, it can be used for the high spin particles and provides 
us with the gauge invariant transition amplitudes. We apply this method to the case 
when the photon is emitted by a constituent in the intermediate state (additive quark 
model). We perform the momentum operator expansion of the spin amplitudes for the 
decay processes. The problem of the nilpotent spin operators is discussed. 

1 Introduction 

In this paper, we provide a consistent presentation of the method of calculation of the 
radiative decays of the composite QQ systems, with arbitrary J . We consider the radiative 
decays which are realised through the transition (QQ)i n — > r y(QQ)out shown in Fig. 1 
(additive quark model approach). The method is based on the spectral integration over 
the masses of composite particles (QQ)i n and {QQ) ou t , it is relativistic and gauge invariant. 
The obtained amplitude is determined by the quark wave functions of the composite systems 
(QQ)in an d (QQ)out- The method, in its substantial part, uses the spin operator expansion 
technique developed previously in [1]. 

The consideration of triangle diagrams in terms of the spectral integral over the mass of 
a composite particle, or interacting system, has a long history. Triangle diagrams appear at 
the rescattering of the three-particle systems, and the energy dependences of corresponding 
amplitudes (on either total energy or one of pair energies) were studied rather long ago, 
though in nonrelativistic approximation, in the dispersion relation technique applied to the 
analysis of the threshold singularities (see [2] and references therein). The relativistic approx- 
imation was used for the extraction of logarithmic singularity of the triangle diagram, e.g., 
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Figure 1: a) Triangle diagram for the radiative transition (QQ)i n — > j(QQ)out with p 2 = M 2 n , 
p' 2 = M 2 ut and (p—p') 2 = q 2 ■ b) Cuttings of the triangle diagram for the double discontinuity 
of the spectral integral with P 2 = s, P' 2 = s' and (P - P') 2 = q 2 . 

see [3, 4] and references therein. Double dispersion relation representation of the triangle 
diagram without accounting for the spin structure was written in [5] . 

In the consideration of radiative decays of the spin particles, one of the most important 
point is a correct construction of the gauge invariant spin operators that allows us to perform 
both the expansion of the decay amplitude (written in terms of external variables) and write 
down the double discontinuity of the spectral integral (written in terms of the composite 
particle constituents). Such a procedure had been realised for the deuteron in [6, 7], corre- 
spondingly, for the elastic scattering and photodisintegration amplitude. A generalization of 
the method for the composite quark systems has been performed in [8, 9, 10]. 

There are two principal points which should be accounted for the processes shown in Fig. 
1 considered in terms of the spectral integration technique: 

(i) The amplitude of the process (QQ)i n — > l{QQ)out should be expanded in a series in 
respect to a full set of spin operators, and this expansion should be done in a uniform way 
for both internal (quark) and external boson states. The spin operators should be orthogonal, 
and the spectral integrals are to be written for the amplitudes related to these orthogonal 
operators. 

(ii) It should be taken into account that in the processes with the real photons (with photon 
four-momentum q 2 — > 0) the nilpotent spin operators appear, their norm being equal to zero 
[11]. Because of that, the representation of the amplitudes may be different, that does not 
affect the calculation result for partial widths. 

The present paper was initiated by the study of the quarkonium systems in terms of 
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the spectral integral Bethe-Salpeter equation — the formulation of this equation for the 
quark-ant iquark systems as well as the discussion of its properties may be found in [12]. 
Since the quark-antiquark interactions at moderately large and large distances cannot be 
considered as well-established, the treatment of quark-antiquark states should be based on 
not only the knowledge on their masses but also on the wave functions, and the source of 
this latter information is the radiative decay processes. Therefore, we present the formulae 
for the radiative transition amplitudes in the form convenient for their simultaneous analysis 
using the spectral integral Bethe-Salpeter equation [12]. 

The paper is organized as follows. 

Section 2 is the introductory one. As an example, we use the transitions QQ(J PC = 
0-+) -> i + QQ(J PC = 1— ) and QQ(J PC = 0++) -> 1 + QQ(J PC = 1~) studied before in 
[8, 9, 10, 11] and recall the method of the amplitude representation in terms of the double 
spectral integral. We also formulate the problem of the momentum expansion of the spin 
amplitude in the case of nilpotent operators. 

In Section 3, the method is applied to the transition QQ(2 ++ ) — > 7 + QQ(1 ) and in 
Section 4 it is applied to QQ(1 ++ ) -> 7 + QQ(1~~). 

In Conclusion, we emphasise that the cases of QQ{2 ++ ) — > , ~f + QQ(l ) and QQ(1 ++ ) — > 
7 + QQ(1 ) are rather general and can be used as a pattern for the consideration of the 
spectral integral representation of the amplitudes (QQ)i n — *• 7 + (QQ)out for the QQ states 
with arbitrary spin. 

2 Radiative transitions P —> j(q)V and S — > j(q)V 

This Section is the introductory one: we consider here the meson radiative transitions 
QQ(0- + ) -> 7(g) + QQ{\—) and QQ(0 ++ ) -> 7(g) + QQ(\—) (below the massive mesons 
with J PC = CT + , 1 ,0 ++ are denoted as P,V,S, correspondingly). These mesons were 
investigated previously within the spectral integration technique [9, 10]. Using these pro- 
cesses as examples, we demonstrate the basic principles of the technique used in subsequent 
sections for the more complicated transitions. In parallel, we discuss the problem of the 
nilpotent operators that arise in the reactions with the photon emission. 

2.1 Transition P -> j(q)V 

Here, we consider the transition P — > 7(g)l / for the virtual photon. We write down the 
spin operator for both initial mesons and quark intermediate states in the triangle diagram, 
with the cuttings shown in Fig. lb. Then, we extract the invariant part of the amplitude 
(form factor) and present it for the emission of real photon expressed through the dispersion 
relation integral. 
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2.1.1 Polarization vectors of the massive vector particle V and photon 
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The polarisations of the vector meson, e^Y\ and virtual photon, e^ 9 ^, are the transverse 



vectors: 

Here, g is the photon four-momentum and p' is that of the vector meson. The pseudoscalar 
meson momentum is denoted as p = q + p'. The polarisation of the vector meson obeys the 
completeness condition as follows: 



9 a p ~ 9af3 -Q- • (/J 



Here, g^p is the metric tensor operating in the space orthogonal to the vector-meson mo- 
mentum p' . 

The polarisation vector of the real photon (q 2 = 0) denoted as e^ 7 -* has two independent 
components only, they are orthogonal to the reaction plane: 

ei 7) g* = 0, cWp a = 0. (3) 
Correspondingly, the completeness condition for the real photon reads: 

- E d 7) W^ )+ W = g£ , 

a=l,2 

£ Qa ' - 9aa' ~ 2 -j" • (4j 

P y± 

Here, g -1 is the orthogonal component of the photon momentum: 

± ± (pi) 

<i a = 9aa'<la' = Qa ~ Pa , 

J_ _ PaPa' / K x 
9aa> - 9aa> ■ \?) 

For virtual photon, (q 2 ^ 0), the completeness condition for polarisation vectors is written 
in three-dimensional space: 

" E #>)# )+ («)=fl#, (6) 

a=l,2,3 

± 7 * _ _ QaQa' 

9aa' ~ 9aa' ^ • 
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2.1.2 Angular momenta of outgoing particles 

The angular momentum of outgoing particles depends on the relative momenta of particles 
in their centre-of-mass system: 



9aa> -TtiV- P')<*' = 9aa> f?"^)^ ?a > ( 7 ) 

Following [1], we determine the angular momenta with the help of operators Xj^).. fiL (q ± ). 
Below, we give an explicit form of these operators for the lowest states considered here: 

L = 0: X(°V) = 1, (8) 

L = 2: X^) = \ { q Ui 2 -\gU<li 

2.1.3 Amplitude for the decay P — > 

The decay amplitude P — > (for the sake of definiteness, the real photon is considered) is 
written as a product of the spin structure and form factor: 

Ap^ 7V - e' a 'ep A a(3 , (9) 

where 

4^ 7V ° = ee a ^ Pu F P ^ v (0) . (10) 

In (10), the electron charge is singled out, and is the wholly antisymmetric tensor. This 
expression works for the virtual photon transition (7^7*) with corresponding substitution: 
Fp^ v (0) -> F P ^v(q 2 ). 

In what follows the important role is played by the spin operator which enters (10): 
^a/3nuQnPu (or, in an abridged form, e a /3 q ± p ). Since e a p vv = 0, the spin operator can be 
represented as follows: 

s { a P ^ v \p,q) = e a p qp . (11) 

Let us stress once more that this spin operator is valid for the reaction with both real and 
virtual photons. 

2.1.4 Partial widths for P -> 7V and V -> 7P 

The partial width for the decay P — > 7V is determined as follows: 

Mpiw = /^ 2 W)|£4^ 7y) l 2 = f^^l^(o)l 2 , (12) 

J a/3 S 1V1 P 
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where My and Mp are the meson masses. The summation is carried out over the photon 
and vector meson polarisations; in the final expression a = e 2 /4n = 1/137. 

The same form factor gives us the partial width for the decay V — > 7P: 

1 f „ , Av^ 7 p) ]2 a M v -M l 

^ al3 1 ~ 24 My 

a/3 V 



MvTy^p = If dMP^,p')\E^r P) \ 2 = S^M^ |Fp ^ ( ° )|2 ■ (13) 



2.1.5 Double spectral integral representation of the triangle diagram 

To derive double spectral integral for the form factor Fp^. 7 y(0), one needs to calculate the 
double discontinuity of the triangle diagram of Fig. lb, where the cuttings are shown by 
dotted lines. In the dispersion representation, the invariant energy in the intermediate state 
differs from those of the initial and final states. Because of that, in the double discontinuity 
P 7^ p and P' 7^ p' . The following requirements are imposed on the momenta in the diagram 
of Fig. lb [6, 8]: 

(jfej + k 2 ) 2 = P 2 > 4m 2 , (k[ + k 2 ) 2 = P' 2 > Am 2 (14) 

at fixed q 2 : 

(P' - P) 2 = (k[ - h) 2 = q 2 . (15) 

Furthermore, in the spirit of the dispersion relation representation, we denote P 2 = s, 
P' 2 = s'. 

When we begin with Feynman diagram, the propagators should be substituted by the 
residues in the poles, that is equivalent to the replacement as follows: (m 2 — k 2 )^ 1 — > 
5(m 2 — k 2 ). Then, the double discontinuity of the amplitude A^^ 1 ^ becomes proportional 
to the three factors: 

disc s disc s < A^ lVi ~ L)) ~ Z P _ lV G P (s)Gy( L )(s) x (16) 

xrf$ 2 (P; h, k 2 )d$ 2 (P'; k[, k' 2 )(2n) 3 2k 20 5 3 (k' 2 - k 2 ) x 



xSp \t l5 {ki + mh^'*CK + m)G ( '' L ' 1 \k'){m - k 2 ) 



The first factor in the right-hand side of (16) includes the vertices: the quark charge factor 
Zp^y (for the one-flavour states Zp^y = cq) as well as transition vertices P — > QQ and 
V — > QQ which are denoted as Gp(s) and G v ^(s') (transition V — > QQ is characterised 
by two angular momenta L = 0,2). 

The second factor includes the space volumes of the two-particle states: d§ 2 (P; ki, k 2 ) 
and d$ 2 (P'; k[, k' 2 ) that correspond to two cuts in the diagram of Fig. lb (the space volume 
is determined in (12)). The factor (2-K) 3 2k 20 5 3 (k' 2 — k 2 ) takes into account the fact that one 
quark line is cut twice. 
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The third factor in (16) is the trace coming from the summation over the quark spin 
states. Since the spin factor in the transition V — * QQ may be of two types (with dominant 
S- or dominant D-wave), we have the following variants for G^' L,J \k'): 

L = 0: Gf^\k') = 7 f\ (17) 
L = 2: Gf^\k') = yfypXfpW . 

Here, k' = (k[ — k 2 )/2 is the momentum of outgoing quarks: k' _L P' — k[ + k 2 . 

The total vertex G^(k') of the vector state is the sum of vertices for the components with 
L = and L = 2: 

G v p (k') = G^WGvwW + G^ l \k')G v(2) {s') (18) 

Correspondingly, there are two traces for two different transitions: P — > 7V(0) and P — > 
7 V(2): 

S P %^ Vm = -5p[G^ 1) (^)(^ 1 + m)^(fci + m)<7fe(-fe + m)] , (19) 
Spg-nrn*)) = _ 5p [GaAi)(^)(^ i + m )^(fc 1 + ro ) i78 (_^ + m )] 

Recall that 7^ 7 * = 7 a'#aa*, see Eq. (6). 

To calculate the invariant form factor Fp^yrn^ 2 ), we should extract from (19) the spin 
factor analogous to S < fg^ lV \q,p) given by Eq. (11). The total form factor is the sum as 
follows: F P ^ v (q 2 ) = F P ^ v{0) (q 2 ) + F P ^ lV{2) (q 2 ). 

For the QQ quark states, this operator reads: 

S { y^\q,P') = e a ^ P ,, (20) 
where q = P' — P, while P' = k[ + k 2 and P = k\ + k 2 . Therefore, we have: 

s P { :r v{L)) = s^^&nsp^vvwrf) , (2i) 



where 



s ~™^'' q) ~ (s%:^-%ns^-%pf (22) 



As a result, we obtain: 

Sp^v(o)(s,s',q 2 ) = Am, (23) 

n i i i\ m \ * n 2 x 6ss'q 2 

S P ^ v(2) (s,s,q) = (2m + s) - 
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with 

A = (s - s') 2 - 2q 2 (s + s') + q\ (24) 
The double discontinuity of the amplitude (16) is equal to 

disc s disc s 'A^ lV{L)) = S^ /3 + ^" yl \q, P') disc s disc s >F P _^ L )(s, s' ', q 2 ) , (25) 

where 

disc s disc s ,F P ^ v{L) = Z P ^ v G P (s)G v{L) (s')d$ 2 (P; h, k 2 )d$ 2 (P'; K, K) x ( 26 ) 

x(2tt) 3 2A; 2 o5 3 (4 - k 2 )S P ^ v{L) (s,s',q 2 ) . 

It defines the form factor through the dispersion integral as follows: 



oo , oo 



, ds f ds' disc s disc s ,F P ^ viL) (s,s',q 2 ) 
W(1)(9 } = J V J - (s - M 2 )(s> - M 2 ) (27) 

4m 2 4m 2 ^ ' 

We have written the expression for Fp^y^^q 2 ) without subtraction terms, assuming that 
the convergence of (27) is guaranteed by the vertices Gp(s) and Gv(l)(s'), see Eq. (26). 
Furthermore, we define the wave functions of the QQ systems: 

/ i \ G P (s) GV( L )(s') 

Ms) = Md(s) = yzrM^ ■ ^ 

After integrating over the momenta, in accordance with (26), one can represent (27) in the 
following form: 

oo 



rp ( 2x 7 [ dsds W{-ssq-m*A{s,s,q' , 2 

Fp^ lV (L){q ) = Zp^ lV[L) J -r^ipp{s)ilj viL) {s ) . S P ^ V (L){s, s , q ) , 

Am? \JK S i S ,q ' 

(29) 

where Q(X) is the step-function: Q(X) = 1 at X > and 0(X) = at X < 0. 
To calculate the integral at small q 2 , we make the substitution: 

s = £ + ^Q, s' = E - ^g, g 2 = -Q 2 , (30) 
thus representing the form factor as follows: 

F^H? 0) = W W / f ^(E)^ (ll (E) / § ^f'-f , 



6 = t /s(4 - 4), A(E, z, Q 2 ) = (z 2 + 4E)Q 2 . (31) 
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After integrating over z and substituting X — > s, the form factors for L = 0, 2 read: 

°f ds s + Js(s — Am 2 ) 

(0) = Z P 

S \ S[ S TiTTi J 



(32) 



4m 2 
oo 



/ds 
—lfjp{s)llJ V ( 2 ){s) X 



X 



4m 2 

(2m 2 + s) In 



Am 2 



3 J s(s — Am 2 ) 



y/s — V s — 4m 2 

The form factors (32) are expressed through ipp(s) and ?/V(L)( s ) by the integral over s: 
these expressions are sufficiently convenient for the combined fitting to the spectral integral 
Bethe-Salpeter equation [12] and radiative decay partial widths. 

The total form factor is the sum as follows: 

Fp^ v (0) = Fp^ v{0) (0) + Fp^ v{2) (0) (33) 

The considered decay P — > offers rather simple case for the representation of the 
process of Fig. 1 in the form of the spectral integral, while the process S — > '-/V gives us a 
more complicated example — there we face the problem of the nilpotent spin operators. 



2.2 Decay of the scalar meson S — > 

In the decay of the ++ meson, the final particles can be in the S- and D-wave states and 
the corresponding spin factors read: 



S-wave : ( e (D e (v> 



D-wave : ■ (34) 



For real photon q^e^ = 0, so, using X^{q L ) from Eq. (8), we see that only the g^p term 
works in the D-wave. As a result, the D-wave operator, 

4V) 4 7) 4 y) = -\ tJ$9&P = -\ , (35) 

gives us the structure similar to that of the ,S-wave. Therefore, the amplitude of the real 
photon emission is determined by a single spin factor in the amplitude: (e^e^). This factor 
may be represented as follows: 

Wg^jP , (36) 

where the metric tensor g^p is given by (4). Thence, we have the following spin operator 
for the amplitude S — > 7V, when the real photon is emitted: 

S%^ v \ P , q ) = g#. (37) 
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This operator takes into consideration both S- and D- waves in the transition S — > 7 V . 

However, the representation of spin operator for the radiative transition S — > 7V is not 
unique: along with (37), one may use a number of other forms. 

2.2.1 Transition amplitude S — > 7V and ambiguities in the representation of 
spin operator 

Using the spin operator (37), we represent the amplitude S — > 7V as follows: 

4fr 7V ° = ^^(p, 9)^5^(0), (38) 

where -Fs^ 7 y(0) is the transition form factor. Below we demonstrate that the representation 
of the decay amplitude (38) is not unique. 

At q 2 = 0, the spin operator (37) may be written as follows: 

p' 2 1 

gip = 9ap(Q) = 9ap + j—^ q a qp ~ (p' a qp + q a p'p) ■ (39) 
Keeping in mind a consistent consideration of the case q 2 7^ 0, here we change the notation 

But it is also possible to apply another spin operator in (38), that is done rather often: 



SSr^ip, q) — Sap = 9aP - H , (40) 



Here, as above, the index a relates to the photon and f3 to vector meson. At q 2 = 0, this 
operator obeys the requirement of transversality: 

q a S afS = , p'pS a p = , (41) 

so it may be equally applied to the S — > 7V amplitude. 

The ambiguity in a representation of the amplitude is due to the existence of the nilpotent 
operator L a p(0), 

L a p(0)L a p(0) = , (42) 

• ^0): 



which is orthogonal to ^^(O): 



9^(0)L a p(0) = . (43) 
The operator L a p(0) obeys the requirement of transversity, 

q*L a p(0) = , L a p(0)p'p = (44) 

This operator can be easily calculated using Eqs. (42) and (43) — it is equal to [11] 

p' 2 1 
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It is seen from (40) and (45) that 

^(0) = S a p + L a p(0) . 

Generally speaking, one can construct the spin operator using any linear combination of 
g^(0) and L a p(0): 

S { a S f jV) = g#(0) + C(p 2 ,p' 2 )L a p(0) . (46) 

Any of these operators may be equally applied to the construction of the transition amplitude 
S — > 'yV for the emission of the real photon. 



2.2.2 Transition S — > in the case of virtual photon q 2 ^ 

For the virtual photon, q 2 ^ 0, the amplitude spin structure is described by the g^p and 
L a p operators which should be generalised for this very case. In terms of p' and q, these 
operators read: 

q 2 

9«0(p',q) = 9a/3 + ^ lq)2 _ p/2q2 PaPp + 

+ iv'J- pw qaq " ~ (P> q y-P>w M + pM ' (47) 



and 



q 2 , , p' 2 



L <*V> q) = {piqY-p'W + {p'qY-p'W g ^ 



^ a rl P ' 2q2 r! a ( AR\ 



(p'q) 2 - p' 2 q 2 ^ [(p'q) 2 - p' 2 q 2 ] (p'q) 

As is easy to see, these operators obey gauge invariance and are orthogonal to each other. 
At q 2 — > 0, the operators g^piP o) an d L a/3 (p', q) transform into the formulae (39) and (45), 
accordingly. 

The transition amplitude S — > j*V is determined by two form factors, namely, 

AafT 1 V) = 9ip{p\q)Ftansverse{q 2 ) + L a p(p ', q) Fi ogitudina i(q 2 ) . (49) 

The first term in (49) corresponds to the process with the transverse polarised final state 
particles (7* and V), while the second one corresponds to the polarisations lying in the 
reaction plane. 

The operators g^pip'o) and L a p(p',q) are singular. To avoid false kinematical singulari- 
ties in the amplitude A^^ 1 v \ the poles in g^p(p'q), L aj3 (p', q) should be cancelled by zeros 
of the amplitude. 
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Instead of Eq. (49), one can use another representation of the decay amplitude, for 
example, that is based on the terms given in (34). Then, we have 

S - wave operator : q^, , 

D - wave operator : q££ X%{q L )g^ . (50) 

Here, g^J is determined by Eq. (3) and gr^J, by (6). Note that the convolution g^, g^ 
does not coincide with g^(p', q) given at q 2 ^ by Eq. (47). There is another, much more 
important difference between the operators given by Eq. (50) and those given by Eqs. (47) 
and (48): the operators (50) are not orthogonal to each other, in contrast to g^{p',q) and 
L a p{p',q). Indeed, the convolution of operators (50) gives: 

& 9%9% X^rtsft = SrT^ + P >2 + <?)■ (51) 

oq p 

Non-orthogonality of the operators (50) is due to the fact that the convolution has been 
performed with the help of the metric tensors g^ an d g~gX ■ Had we operated with the 
normal metric tensor, namely, had we substituted in (51) g^2> — > g aa i and g^Y, gppi, 
we would have the orthogonal S- and .D-wave operators. The metric tensors g^2> and 
g-pp in (50) allow us to fulfill the gauge invariance — in this way, just due to the gauge 
invariance, the orthogonality in the S- and D-w&ve operators (50) is broken. But, in the 
spectral representation of form factors of the composite systems, the orthogonal operators 
are needed to avoid the double counting. This is the reason why furthermore we deal with 
the orthogonal operators represented by formulae (47) and (48). 



2.2.3 Analytical properties of the amplitude at q 2 = 

Let us turn to the real photon and discuss analytical properties of the amplitudes, namely, 
the cancellation of kinematical singularities. The A^p*" 1 amplitude reads: 



9af3 + 



Am 2 , 



mt 



rriy) 2 



QaQ/3 



m 



v 



iPaQ/3 + QaPp) 



transverse 



(0) + 



+ 



4m? 



(m| - m\Y qaQp 



o 9 Q&Pli 

rrig — my ' 



longitudinal 



(0) 



(52) 



Here, we have used 2{p'q) = m| — m v . To make nonsingular the term in front of q a qp at 
m| — > my, it is necessary that 



[Ft 



(0) + F longitudina i(0)] m s i 2 ~ (m 2 s -my) 2 



transverse 



(53) 



This requirement is sufficient for the cancellation of kinematical singularity in front of q a p'p- 
However, to remove kinematical singularity in the term p' a qp, the following condition for 
F t ransverse{fy should be fulfilled: 



transverse 



(0) 



(m 2 s — my) at (m 2 s — m v ) — > . 



(54) 
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The constraint (53) is in fact the requirement imposed on F longitudinal (0), but the F longitudinal (0) 
itself, as was noted above, does not participate in the definition of the decay partial width 
of the process S — > 7V. 

The second constraint given by Eq. (54) for F transverse (0) is the principal one for the 
decay physics — in quantum mechanics it is known as Siegert's theorem [13], see also the 
discussion in [14, 15, 16]. 



2.2.4 Spin operator decomposition of the quark states in the triangle diagram 

Now, let us consider the quark triangle diagram of Fig. 1. As was said above, in the QQ 
systems there are two possibilities to construct vector mesons with the angular momenta 
L = and L = 2. For the transitions V — > QQ(L), we apply the vertices introduced in (17): 
G^' ' 1 ^ and G^' 2 ' 1 ^. For the transition S — > QQ(L), we use, following [1], the spin operator 
ml, where / is the unit matrix. The traces for two processes with the different vector-meson 
wave functions (L = 0,2) read: 



Sp K ap W = ~M G f3 (K + + m)ml{-k 2 + m)\ , (55) 

: , s r vm = -s P [Gf 



Spl S p V(2)) = -Sp[G { } A1) (ki + m)^*(ki + m)ml(-k 2 + m)] 



To calculate the invariant form factor for the transverse polarised final state particles (we de- 
note it in an abridged form as Fs^y^^q 2 )), one should extract from (55) the corresponding 
spin factor. For the quark states, this operator reads: 

4 ++ ^ 7l " ) (^n = g^n- (56) 

Recall that P = k[ + k 2 and q = P - P' = £4 - k[ . We have: 

SpTr ViL)) = S% ++ ^--\q,P')S s ^v {L) (s, S ',q>), (57) 



where 



The spin factors Ss->-yv(L)( s , s', q 2 ) at L = 0, 2 are equal to 

Hs,s',q 2 Y ' 



S s ^v(o)(s, s', q 2 ) = -2m[(s - s' + q 2 + 4m 2 ) - / S f ] , (59) 



Ss^v(2){s,s',q 2 ) = -^- 7 =[4m 4 -2m 2 (3s + s , -g 2 ) + s(s-s / + g 2 ) + 
+ w 2SSV 2 , (16m 2 + 3g 2 - S -3/)], 
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with A(s, s', q 2 ) given by (24). 

The form factor of the considered process takes the form: 



oo 



jp i 2n 7 f dsds' ,6(-ss'q 2 -m 2 \(s,s\q 2 )) , 2 

F s ^v(L){q ) =Z s ^v J -r^^s{s)ipv(L){s) , = S s ^ lV (L)(s, s , q ) . 

4m 2 yA(s,s ,g ) 

(60) 

To calculate the integral at q 2 — > 0, we make, in a complete similarity with the calculations 
done in Eqs. (30) and (31), the following substitution: q 2 = —Q 2 , s = S + zQ/2, s' = 
£ — 2<5/2 and represent the form factor at small Q as follows: 



oo +6 



F S .W-Q 2 - 0) = W / f «£)^>(£) / f ^fe^ , (61) 

4m 2 -6 16^A(S,Z,g 2 ) 



where b = y £(£/m 2 — 4) and A(£, z, Q 2 ) = (z 2 + 4S)Q 2 . After the integration over z and 
substitution E — > s, we have: 

s)^v(o)(s)I s ^v(s), (62) 





oo 

m r 




4^ y 


7T 




4m 2 






oo 






m y 




27 J 


7T 



4m 2 



/ x n ; — 7T ^ 9 , \fs + Vs — 4m 2 

Is^v{s) = Js(s - Am 2 ) - 2m 2 In v V= r=? ' 



The total form factor is equal to 

(0) = F S (0) + F S (0) . (63) 

2.2.5 Partial widths for the decay processes with the emission of real photon 

Similarly to the form factor calculations performed above, the partial width of the scalar 
meson decay S — > reads: 



MsTs^v = /*Mp;?,jOIE^°I 2 = f |F^(0)| 2 . (64) 

a/3 S 



Recall that in the final expression a = e 2 /An = 1/137. Likewise, partial width of the vector 
meson decay V — > is equal to: 

MvTy^s = ^IW«)| 2 (65) 
Here, we do not specify the quark structure of the vector meson omitting the index L. 
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2.2.6 Normalization conditions for wave functions of QQ states 

It is convenient to write the normalization conditions for P, S and V meson wave functions 
using the charge form factor of a meson: 



charge 



(0) = 1 . (66) 



The amplitude of the charge factor is defined by the diagram of Fig. la, with (QQ)i n = 
(QQ)out- For P and S mesons, the amplitude takes the form: 

Ax(q) = e(p + p) a F charge (q 2 ) , (67) 

while F c h arge (q 2 ) can be calculated in the same way as the transition form factors considered 
above. For the meson V, we consider the amplitude averaged over the spins of the massive 
vector particle. At q 2 = 0, it takes a form: 

A%l(q - 0) = 3e(p + p') a F^ rge (0) . (68) 
The normalisation conditions based on the formula (67) for P and S mesons read: 



1 - /i^W 1 ^ < 69 » 

4m 2 
oo 



1 = / i^M 2m Ms- 4 ™ 2 )\F^- 

4m 2 

For the vector mesons V, the normalisation condition reads as follows: 

1 = W 00 [V] + W 02 [V] + W 22 [V], (70) 

WnW] = I J ^ ^(o)(*) 4 (* + 2 - 2 ) \J S -^ . (71) 

4m 2 

W 02 [V] = ^ J ^ Mo)(s)M2)(s) (s - 4m 2 ) 2 J S -^- , 



4m 2 



2 °f ds 2 (8m 2 + s)(s - 4m 2 ) 2 / s - 4m 2 
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4m 2 

In [9, 10], the calculations of charge form factors are explained in a more detail. 



3 Transitions 2 ++ -meson —> and l ++ -meson —> 

In this Section, using the decays of the 2 ++ -meson (denoted as T) and l ++ -meson (denoted 
as A), we perform the treatment which can be easily generalised for any spin particle. 
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3.1 Transition T — > 

To operate with the 2 + -meson, we use the polarisation tensor e^ u (a) with five components a = 
1, . . . , 5. This polarisation tensor, being symmetrical and traceless, obeys the completeness 
condition: 

1 / 2 \ / / 

E e F( fl )^v( a ) = 2 [dnn'Qvu 1 + 9[jlu'9uh' ~ ~^9nu9n>v>) = , 

E ^(a)e+ (a) = 5 . (72) 

a=l,...,5 

Here, Oj^ is a standard projection operator 0^ v , for the angular momentum L = 2 which 
obeys the requirements: O^^nO^,", = 0^ v , and O^, = 0, see [1] for more detail. 

In terms of the polarisation tensor and vectors e£*\ the independent spin struc- 
tures for the transitions with virtual photon (q 2 ^ 0) read: 

(1) S- wave : t^e^^p , 

(2) D-wave: ^(g^e^''^) , 

(3) D-wave: e^g^^ef ^ , 

(4) ^-wave: rffe 1 )^'^ , 

(5) G-wave: e^^^e^ . (73) 

Correspondingly, we have five independent form factors which describe the transition 2 ++ - 
meson — > 7*V. But, for the real photon (g 2 = 0), the number of independent form factors 
is reduced to three and, keeping this fact in mind, we consider below the production of the 
transverse polarised photon. 

3.1.1 Spin operators for transverse polarised photon, 2 + -meson — > 7* ± V 

Here, we consider the transverse polarised photon, though with q 2 7^ 0. We introduce the 
following spin operators, which correspond to the spin structures (73): 

(74) 



c(l) 
fiv,af3 




o(2) 

^ [XV,afi 


'1 


o(3) 
ixv,af3 


1 fy v X r(2) fn-hn- L - L n- Ly 
9 °V ^V/J'W )9 l x'a9f}>(j i 


o(4) 
^ fiv,af3 


1 (l^'^X^ frr L ^- LL n- Ly 


o(5) 


IncVyW / ±\ ±± ±v 

4 ^1/ A /iVa'|3'W )9 a 'a9/3'(3 ■ 
Q± 
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Recall that here q£ = g^ a ,q a > = q a - p a (pq)/p 2 and g^ a , = g aa , - p a p a >/p 2 . 

Let us demonstrate the method of construction of these operators by considering the 
G-wave spin structure from (73), V^'^u^'a'/3'(?" L ) e cJ ^ e %' '■ ^ should be multiplied by the 
polarisations e+ u (a), e^* ^ (b) , e^ + (c), with the summation over a, b, c: 

E^(«)v^(«)^^(? ± )^ ± H6)cr ±)+ (6)^( c )4 m ( c ). 

a,b,c 

that gives us S^j a/3 , because the convolution e^> ± \b)e^ / "'- L " >+ (b) results in g^. 
The operators (74) should be orthogonalised as follows: 



o(±3) / / \ o(3) o(±l) / / \ (^ vU'b^^lv.a'/j') 

^ ^u,a/3\P y 1 ) — ^Liv,af3 ^fMU,a/3\P ) / (±1) 



In this way, we construct three operators, i = 1,2,3. The operators Sj^^p an d ^it,a/3 are 
nilpotent at g 2 = 0, so we do not present explicit expressions for them here but concentrate 
our attention on the calculation of the amplitude for the emission of the real photon. 

The results of calculation of the spin operator convolutions are presented in Appendix 2. 
Here, we give the orthogonalised operator norms, which determine the decay partial width: 

S^%{p\q)S^%{p\q) = zi(MlM*,q 2 ), (76) 

, 3M| + 34M|M 2 + 3M V 
z u (M T ,M v ,0) = , 

7 ± (M 2 M 2 m Mj + 10M|M y + My 

Z 22 (M T ,7W y ,Uj - y 3M 4 +34M 2 M 2 +3M 4 ' 

zHU 2 M 2 0) - 9 W + ^) 2 

3.1.2 Calculation of the transition amplitude T(L) — > 7V(I/) for the emission of 
the real photon 

The transition amplitude of the T — > 7V decay can be written using the operators (75) as 
follows: 

i=l,2,3 
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where Ft(l)^v{u)(®) are the form factors at q 2 = 0. 

Considering the double discontinuity related to Fig. lb, we should expand the following 
traces in a series over the spin operators: 



(T(l)-nrV(O)) 



q^ 1 ^ 



o (T(3)-KyV(0)) 
^ .r fiv, a/3 

o (T(3)-KyV(2)) 



----- -5p [^ 1A1) (^ + m) 7 ^(fci + m)T^(k)(-k 2 + m 

(k\ + m)^(k\ + m)TW(fc)(-4 + m 

,(3) 

/it" 



G 



(1,2,1) 



(78) 



-Sp [G^ 0,1) (fci + m) 7 ^(fci + m)T^(k)(-k 2 + m) 
-Sp [G^ 2 ' l \k\ + m)7a 7 (£i + m)T$(k)(-k 2 + m) 



Recall that here we have used the notations 7^ 7 = g^la'ik = (k 1 —k 2 )/2, k' = (k[ — k 2 )/2 
and vertices df^\ df 2 ' l) given in Eq. (17). 

The operators for the transitions 2 +4 -meson — > QQ(L) for L = 1, 3 read: 



5 

71 



2 

^7^ + &i/7/i ~~ 



(79) 



To calculate invariant form factors F^l_^ v , L Jq 2 ), we should expand the traces of (78) in a 

series over the spin operators for the quark states S^l j3 (P',q) and extract the invariant 
factors: 



i=l,2,3 



(80) 



Let us emphasise that in (80) the spin operators depend on the intermediate-state quark 
variables, P' and q. 

The invariant spin factors read: 



{ S P. 



(T(L)^V(L')) s (±i) 



P', q 



{P'A)S^UP'A) 



(81) 



where i = 1,2,3. Invariant spin factors determine the form factors in a standard way: 



^r(L)-^ T y(L')(^ 2 ) 



4m 2 



dsds' 
16tt 2 



^T(L)(s)^y(L')( S ') X 



(82) 



g(-ssV-m 2 A(s, s / ,g 2 )) o(±t) , 2 

x = °r(L)^ 7 y(L')^ ' ' 



v/A(s,s', g 2 ) 

To calculate the integral at g 2 — > 0, we make, as before (see Eqs. (30) and (31)), the following 
substitution: q 2 = —Q 2 , s = £ + zQ/2, s' = £ — zQ/2 and perform the integration over z. 
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We have: 



Here, 



where 



F (i) 



= z 



? (1) , 



>T(1)^V(0) 



(s 



S, 



s, 



'T(1)-»7V(2) 

? (2) , 
T(1)-»7V(2)* 

? (1) ( 
>T(3)-*7V(0) » 

(2) , 
T(3)-KyV(0)l 

>T(3)-KyV(2) 

(2) , 
T(3)^ 7 V(2)l 



' L ) -,H Z- ' ) / ^2 ^ T(L) (i0 ^ J T(L)^Vi L ■ 

Am 2 



^(8m 2 + 3 S )4V( s ) . 



Ho(3) 

3°T(1)- 7 V(0) 

40 



-9 (3) M 



3^ 



r(2) 



_^_(l 6m 2 _ 3s)(4m 2 _ s )/(^ v ( a ) , 



V2 



20 



3 D T(3)^7V(0)l 6 



► 7 VV 

v/2, 



r(2) 



(3) 



3 ^T(3)^ 7 V(2) 



» 



1 

"72 



(16m 2 - 3s)(4m 2 - s)l£> v (s) , 



l£l jV (s) = 2m 2 ln^±^g-,/^^), 



r (2) 

-IT 



m 2 (m 2 + s) In 



y^s + Vs — 4m 2 1 



Total form factor is a sum over four terms: 



^s(s-4m 2 )(s + 26m 2 ) 



2^ r T(L)- 7 V(L')' 



L,L' 



3.1.3 Normalisation conditions and partial widths 

The normalisation condition for tensor mesons reads: 
1 = W 11 [T]+W 13 [T]+W 33 [T], 
W U [T] ' > " -'- 2 '-^'o-^ 



; /• | / 4m 2 

5 7 T ^^( 1) («)2(8m 2 + 3 S )( S -4m 2 )^l-^ 

4m 2 
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W 13 [T] -- 

pi - 17 



1 °° d \/3 

4m 2 



1- 



4m 2 



4m 2 



Partial width of the decay T — > 7V is equal to: 

/iu,a/3 

a m\ — my 



20 m| 



4(m|,m^,o) (4!V(o)) 2 + 

+ 4(M|,M y ,0) (4 2 _i 7y (0)) 2 + 4(M|,M2,0) (41^(0))' 
The same block of form factors determines the partial width for V — > 7T: 

/if, a/3 

4(714,714, 0) (4V(0)) 2 + 
+ 4(i4,M y ,0) (4 2 i 7F (0)) 2 + 4(M2,M2,0) (41 7F (0)) ; 



m y r y ^ 7r = e 2 J d<& 2 {p],q,p')\ ) I^W/sl 



a m 2 / — to 2 - 
12 to^ 



(89) 



(90) 



Let us emphasize that the factors z^ a (M^, M v ,0) are symmetrical with respect to T <-> V 
permutation: 4 a (M 2 ,M y ,0) = 4 a (M y ,M 2 ,0). 



3.2 Transition A -> 7"!/ 

For the reaction l ++ -meson — > 7I -meson (or A — > 7V), one can write three partial 
states: the S'-wave state with the total spin S = s 7 + Sy = 1 and two .D-wave states with 
S = s 7 + sy = 1 . Generally, we have three spin structures, but only two of them survive in 
the case of transverse polarised photon (below, as before, p is the momentum of the decaying 
particle and q is that of the outgoing photon): 



cW 1 \ _ ±± ±v 

^UtapyPiQ) — 9 aa' 9 1313' £ ^a'/3'p , 



(91) 



9 1f3'9^'9a^9pp' Zfi'a'q^p 



L _L ±V 

' 9 C Lf3'9f3(3' £ fiaq- L p , 



S ^ap(Pi Q) - ~^Qa'9w'9aa'9pi3' £ ^'(3'q ± p ~ ■ 

Q± 

(3) 

Here, as previously, we use the abridged form e^ a ^p^ = e m p v . The vanishing of S^ aj3 is 
due to the equality qfg^ = 0. 
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3.2.1 Spin operators and decay amplitude 

The operators S^ a/3 (p, q) should be orthogonalised: 



S^( P ,q) = S^(p,q), (92) 



,(±2), _ C (2) , s C (_:l.) ^ i^li',a'P'^P'> Q^^'li'la'P'^Pi 



We determine the convolutions 

0> ?) = ^(^1, My, g 2 ) . (93) 

At q 2 = (see Appendix 3 for the details), they are equal to: 



ztl(M > M = zMlWtM , (94) 



(Mj + 6MiMy + My) 

The transition amplitude A — > 7V reads: 



4£y v(£,)) = E^(f-^U(°)- (95) 



i=l,2 



being determined by two form factors v (0) (i = 1,2). 



3.2.2 Calculation of the quark triangle diagram of Fig. 1 for the emission of 
the real photon 

The diagram of Fig. lb for the processes A — > 'yV(L) (L = 0, 2) is determined by the 
following traces: 



o (A-^V(0)) 

where the transition A 



—Sp 



G { ^' l \k\ + m)^(h + m)T^k)(-k 2 + m) 



(96) 



: -Sp [Gf 2 ' l \k\ + m)^(k\ + m)T^(k)(-k 2 + to) 
QQ is equal to: 

'2 



T,{k) 



(97) 



Recall that k — (ki — k 2 )/2 and P — k± + k 2 . 

To calculate the invariant form factor F4^ 7 y( L )(g 2 ), we should expand (96) in a series 
with respect to the spin operators Sjj^^P, q) where q = P — P'\ 

SP ( ^ V{L)) = ZsWftfis^y^sW), (98) 

i=l,2 
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that gives us 

where i = 1 or 2. 

The transition form factor is determined by the standard formula: 

^itW* 2 ) = / ^a(^(L)M x (100) 

4m 2 

e(-ss'q 2 - m 2 X(s,s',q 2 )) ^ , 



In the limit q 2 — > 0, performing calculations as in the cases considered above, we obtain: 



^ilrVW = Z a^v J ^(8)M*)ML)(*)tfLrV{L)(*) . ( 101 ) 



'3 



4 2 i 7 y(2)( s ) = — (4m 2 -s)/^ 7 y(s), 



where 



W^) = f2m 2 In ^ - " 4m 2 )) . (102) 

V v s ~~ V s — 4m/ / 

The total form factor is equal to: 

(0) = F A (0) + F A _ >7 y(2)(0) (103) 

3.2.3 Normalisation condition and partial widths 

The normalisation condition for the 1 ++ meson wave function reads: 



^/lil^M^-W 1 -^- (104) 



4m 2 

The partial width of the decay A — > 7V is equal to 



m^r^y = e 2 J d$ 2 {p; ,q,p')^ IA»,a/sl 2 ( 105 ) 

4(M 2 ,M 2 ,0) (FW(0)) 2 + 4(M^,J^,0) (f^o)) 1 
22 



a m 2 ^ — m 2 / 
12 ml 



For the partial width of the decay V — > ^A, one has: 




,-L 
11 



(M v ,M 2 A ,0) (F( 1 )(0)) 2 + 



(106) 



12 m?, 



+ 4(M^,Ml,0) (^(O)) 2 



Let us emphasise that ^(My, M\, 0) ^ ^ a (M|, My, 0). 

4 Conclusion 

The considerations of the tensor meson decays (Section 3) and pseudovector meson decay 
(Section 4), being in fact general cases, can be easily expanded for any mesons. 

Actually, the tensor meson decay is a pattern for an amplitude, where the parity of initial 
meson coincides with the parity of final state. For this case, we construct the spin scalars 
from the polarisation and angular momentum functions Xf£). L (k ± ), see Eq. (73) for tensor 
meson. With a completness condition for the vector and tensor polarisations, we construct 
gauge invariant spin operators ((74) for the tensor mesons. The orthogonalisation of these 
operators for the case of the real photon emission allows us to single out the operators with 
nonzero norm, of the type of Eq. (75)), and nilpotent operators. The operators of the first 
kind are used in the expansion of the amplitude in a series in respect to external operators 
(Eq. (77)), as well as for the quark triangle diagram (Eqs. (80) and (81)). The spectral 
integrals are written for the invariant form factors, which are the coefficients in front of the 
orthogonalised operators. 

The decay of axial meson provides us with an example, where the parity changes sign 
from initial to final state. In this case, there is only one difference: in the first step, we 
construct the pseudoscalars from the particle polarisations and angular momentum function 
Xf£) (k- 1 ). Further consideration is carried out following the scheme common to that of 
tensor mesons. 
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Appendix 1: Convolutions of spin operators 
in the decay T —> jV 



S%, a0 (P', Q)S { Xs(P'. 9) = s', q 2 ) 



■Q' 



■ q' 



In the calculation of the spin factors we need the convolutions 

They are as follows: 

2n (s, s 

Z22(S, S 

z u (s, S 
z 12 (s,s 
z 13 (s, s 

z u (s, s 

Z2A (S, S 

z u (s, S 

215 (S, S 
^25( s 5 s 
z 35 ( s , s 



■ q' 



■ q' 



■ q' 



■ q' 



■ q' 



z 45 (s, s',q 



For the orthogonal operators, we have: 





;i.s- 


- 6sq 2 + 34ss' + 3g 4 - 6q 2 s' + 3s' 2 


= 




1 2ss' 




d J 






o„2 
OS 


^^^,2 i iQ.J i q~4 £?~2^/ i q ^/2 

— osf/ + loss +05 — 05 S + OS 






izss 




OS 


c^>2 _j_ Q/| „„/ _j_ „4 „2 „/ i „/2 

— usq + o4ss + oq — vq s + os 






48ss 




1 

1 , 






o 

OS 


2 c 2 / i 4 c 2 / i /2 

— osg — 8ss + 3g — og s + 3s 


= 




12ss' 




3s 2 


- 6sg 2 + 34ss' + 3g 4 - 6g 2 s' + 3s' 2 


— 


1 


24ss' 




2 ' 






i 
i 






2 






-3s 


2 + 6sg 2 + 8ss' - 3g 4 + 6g 2 s' - 3s' 2 






24ss' 




2s 2 


- 4sg 2 + Wss' + 2g 4 - 4g 2 s' + 2s' 2 






8ss' 




-s 2 


+ 2sg 2 + 2ss' - q 4 + 2g 2 s' - s' 2 




3 


4ss' 




2 ' 






s 2 - 


2sg 2 + 4ss' + q 4 - 2g 2 s' + s' 2 






4ss' 




s 2 


- 2sg 2 - 2ss' + q 4 - 2q 2 s' + s' 2 






8ss' 



S { ^%{P', q)SH%(P', q) = z^ b (s, s', q 2 ) 



with 



zii{s,s',q 2 ) 



3s 2 - 6sg 2 + 34ss' + 3g 4 - 6g 2 s' + 3s 



12 



\2ss' 
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, , 2x s 2 - 2sq 2 + 1(W + q 4 - 2q 2 s' + s' 2 

z 2 2{s,s,q) = 9; 



3s 2 - 6sq 2 + 3Ass' + 3g 4 - 6gV + 3s' 2 



± , 2 9 ( S -g 2 + S ') 2 



2 s 2 - 2sg 2 + lOss' + q 4 - 2q 2 s> + s' 2 
Spin factors given by S, z variables in the limit Q — > read: 

,(i) ^ ~ ,n .ha 3 v / 2(4m 2 S + mV-S 2 )(8m 2 + 3S) 



^W(0)( S ^^-°) = ( 112 ) 



5 T( ) i)-H y v(o)( i: ' 2; ' < 3- > °) - 3 (4S + ^ 2 ) 2 

(1) . , 3 (4m 2 S + mV - E 2 )(16m 2 - 3E)(4m 2 - S) 
^(1)^(2)1^^' v -> UJ - 2 q 4S + z 2 

(2) 2 (4m 2 S + mV - T 2 + S^ 2 )(4m 2 S + mV - S 2 )(8m 2 + S) 
b T{i)^v{2) (^>*> V -»■ U J - 3 (4S + z 2 ) 2 ' 

(1) 3y^ (4m 2 S + mV-S 2 )(4m 2 -S) 2 
^r(3)-» 7 v(o)^> z > V -»■ "J ~ io 4S + 2 2 

(2) 4^2 (4m 2 S + mV - T 2 + E,z 2 )(4m 2 E + mV - S 2 )(6m 2 + S) 
b T(3) ^ v{Q) {h,z,Q ^0) - — (4S + z 2 ) 2 ' 

(1) 3 (4m 2 S + mV-S 2 )(8m ; 

q( 2 ) ,y r n _> m (4m 2 S + mV-S 2 + S^ 2 )(4;.. 

^(3)^(2)1^^ -> uj - 9(4S + ^ 2 ) 2 

(3) 2 (2) 

Appendix 2: Convolutions of spin operators 
in the decay A — > 7V 

We denote the convolutions of the spin operators at q 2 ^ as follows: 

Q)tfU P ', 9) = s '> ^) > ( 113 ) 

where 

, 2A -(s 2 + 6ss' + s /2 ) + g 2 (2s + 2s'-g 2 ) 



8^ (4m 2 S + m 2 z 2 - S 2 + S^ 2 )(4m 2 S + mV - S 2 ) 



(1) 3 (4m 2 £ + mV-£ 2 )(8m 2 + £)(4m 2 -£) 2 
)-n-y(2)(^' z -> Q ~^ 0) — ' 

(2 ) ^ ^ n , (4m 2 S + mV - S 2 + S^ 2 )(4m 2 S + m 2 z 2 - S 2 )(12m 2 - E)(E - 4m 2 ^ 



2s' 



2n(s,s',5) = 

Z22{S,S,q) = — , 

zi2(s,s',q 2 ) = {S + S ^ q2)2 ■ (H4) 
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For the orthogonalised operators, 
one obtains: 

-(s 2 + 6ss' + s' 2 ) + q 2 {2s + 2s' - q 2 ) 



(115) 



^n(s,s',g 2 ) = 

*22(*» s '»0 2 ) = 



2s' 

2s(s + s' - g 2 ) 2 



^s 2 + 6ss> + s' 2 ) - q 2 (2s + 2s' - q 2 ) ' 



(116) 



For the calculation of the spin factors in the double dispersion integral, we use variables 
s = E + \zQ, s' = E — \zQ, and q 2 = —Q 2 . Using these variables, in the limit Q — > 0, we 
have for spin factors: 



.4^7^(0) 



s 



(2) 

A-> 7 y(o) 



o(2) 

°A-> 7 y(2) 
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x E(z 2 + 4E) 

1 (32m 4 E + 8m V + 4m 2 E 2 + 7m 2 E^ 2 - 3E 3 ) 



2^/6 



V^(z 2 + 4E) 



. 
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